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I Abstract Correlation functions of the two-dimensional Ising model on the 

'■^ ■ periodic lattice can be expressed in terms of form factors — matrix elements 

' of the spin operator in the basis of common eigenstates of the transfer ma- 

trix and translation operator. Free-fermion structure of the model implies 
' (— l'' that any multiparticle form factor is given by the pfaffian of a matrix con- 

^ ^. structed from the two-particle ones. Crossed two-particle form factors can be 

obtained by inverting a block of the matrix of linear transformation induced 
on fermions by the spin conjugation. We show that the corresponding matrix 
. is of elliptic Cauchy type and use this observation to solve the inversion prob- 

lem explicitly. Non-crossed two-particle form factors are then obtained using 
theta functional interpolation formulas. This gives a new simple proof of the 
factorized formulas for periodic Ising form factors, conjectured by A. Bugrij 
I and one of the authors. 
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1 Introduction 

The two-dimensional Ising model is the foundation of our understanding 
of critical phenomena and quantum field theory. Starting from Onsager's 
calculation of the partition function [23 , its study has led to remarkable 
developments in many areas of mathematical physics, including the theory 
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of Toeplitz matrices, quantum integrable systems, Painleve equations, and 
conformal field theory. In an effort to keep the bibliography of reasonable size, 
we refer the reader to P^[201[^ for an account of some of these developments 
and further references. 

Correlation functions of the Ising model on an M x square lattice 
with periodic or antiperiodic boundary conditions in each direction can be 
expressed in terms of elementary blocks, called form factors. These blocks are 
matrix elements of the spin operator in the basis of common eigenstates of the 
2^ X 2^ transfer matrix and the operator of translations. Compact factorized 
formulas for the form factors of Ising spin have been conjectured in [SJ[7] on 
the basis of (i) long-distance expansions of the two-point correlation function 
on the cylinder {M — )■ oo) [S] and (ii) direct transfer matrix calculations for 
small values of N. 

One can try to prove the conjectures of [BJ[7| in several ways, reflecting 
different facets of the integrable structure of the Ising model. First, it can 
be seen as a special case of the eight-vertex model. This is not very helpful, 
as the methods of computation of finite-lattice form factors for models with 
elliptic i?-matrix have not yet been developed. Another option is to con- 
sider the Ising model as a particular case of the Baxter-Bazhanov-Stroganov 
T^^^ -model [T1I51IT7]. The corresponding i?- matrix is trigonometric, but the 
i-operator intertwines spin-^ and cyclic type evaluation representations of 

the quantum afhne algebra Uq{sl2)- For this reason, the transfer matrix can 
not be diagonalized by the standard algebraic Bethe ansatz technique, the 
corresponding quantum inverse scattering problem does not admit a simple 
solution, and the Lyon method of computation of form factors for XXZ-type 
models 16^ can not be applied. 

These difficulties can be partially overcome in the framework of Sklyanin's 
method of separation of variables |30j . Local spin operators are then ex- 
pressed in terms of the quantum monodromy matrix by means of an iterative 
procedure. The transfer matrix eigenstates are represented by linear combi- 
nations of eigenvectors of an auxiliary problem, and the computation of form 
factors involves summations over this auxiliary basis. Even in the Ising case, 
performing the sums explicitly and showing that they do reduce to factorized 
expressions of [5J[3 is a highly nontrivial problem, which was solved in [MTU] . 
In spite of this success, the whole procedure looks overwhelmingly cumber- 
some as compared to final answer, and there remains a strong feeling that 
the result can (and has to) be obtained in a much simpler way. 

The most natural approach would be to take advantage of the free-fermion 
structure of the Ising model. This structure was discovered in |15| . where it 
was used to give an alternative derivation of the partition function. 

Recall that the periodic Ising transfer matrix V commutes with the global 
spin flip operator U . Since f/^ = 1, the eigenstates of V split into two sets, 
corresponding to -1-1 and —1 eigenvalues of U . The transfer matrix can be 
naturally written as 




(1.1) 
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where Va^p commute with the projectors 



1±U 
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and are expressed in terms 



of the generators of a CHfford algebra. They actually belong to the Clifford 
group, i.e. the generators transform linearly under conjugation by Va.p- We 
emphasize that the same is not true for the full transfer matrix V. 

Clifford group elements can, in principle, be reconstructed from the in- 
duced linear transformations, see e.g. Chapter II in [3]. The case of Va^p 
is especially simple, as these matrices act on the generators as products of 
commuting two-dimensional rotations; this becomes manifest after discrete 
Fourier transforms (antiperiodic for Va and periodic for Vp). Diagonaliza- 
tion of these rotations gives the spectrum of 14, p and hence, by that 
of V. The eigenvectors of Va and Vp are given by multiparticle fermionic 
Fock states, constructed using two different sets of the creation-annihilation 
operators ^15i . 

Ising spin operator a also belongs to the Clifford group. However, it 
changes the value of Z2-charge U from to —1 and vice versa. Accord- 
ing to (II. ip . these values are associated with fermionic eigenstates of V of 
different types, and it is not clear how one can compute form factors between 
them. The problem can be circumvented on the infinite lattice (N oo), 
where the distinction between two types of fermions effectively disappears. 
In that case, two-particle Ising spin form factors have been found in j29) . 
Theorems 5.2.1-5.2.3, and multiparticle ones in 27], Theorem 5.0. 

Further progress for finite N was made by Hystad and Palmer in two 
recent papers |11II26) . It was noticed that if we combine the creation-annihi- 
lation operators of a- and p- fermions into A'^-dimensional column vectors ip ^ , 
ip, cp'^, (f, and express the result of conjugation of p-fermions by a in the 
a-basis. 



then form factors of a between Fock states of different types can be expressed 
in terms of induced rotations in essentially the same way as if we had an 
ordinary Bogoliubov transformation involving fermions of only one type. For 
example, for T <Tc 

— an arbitrary n-particle form factor is given by the pfafSan of an n x n 
matrix, constructed from the two-particle ones; 

— normalized two-particle form factors coincide with the elements oi N x N 
matrices D'^, BD'^, D-^C. 

Explicit form of the blocks A, B, C, D can be fixed in a rather straightforward 
way. The computation of finite-lattice form factors therefore reduces to the 
inversion of D and calculation of the matrix products BD^^ and D^^C . The 
present work is devoted to the solution of these two problems. 

The key point of our analysis is an elliptic parametrization of the Ising 
spectral curve, equivalent to the elliptic substitutions used in Yang's deriva- 
tion of the spontaneous magnetization [33]. We will show that in this pa- 
rametrization D is given by an elliptic Cauchy matrix. Its determinant and 
inverse can be computed using the so-called Frobenius determinant formula 
[8]. (Let us mention that the Frobenius determinant and its generalizations 
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have recently appeared in the computation of the partition function of the 
eight- vertex SOS model with domain wall boundary conditions PHI^ ). In 
fact taking the iV — )■ c» limit of det D provides a useful alternative to Yang's 
derivation. Further, the products BD^^, D~^C can be calculated using theta 
functional analogs of the Lagrange interpolation identities. The matrix of 
two-particle form factors also turns out to be of elliptic Cauchy type. Com- 
puting its pfaffian and going back to the usual trigonometric parametrization, 
we obtain the general formula for finite-lattice form factors of Ising spin, con- 
jectured in [61I2- 

It is worth noticing that in a recent paper |12j form factors of the quantum 
Ising chain in a transverse field were derived. They represent a limiting case 
of the form factors obtained in the present work. Although the derivation in 
[H] also uses the fermion algebra, it is motivated by the calculation of form 
factors for superintegrable chiral Potts quantum chain [T3| using Baxter's 
extension [2] of Onsager algebra and it does not use linear transformations 
induced by the spin operator. 

This paper is organized as follows. In Section [21 we introduce basic no- 
tation and recall Kaufman's fermionic approach [TS] to the diagonalization 
of the Ising transfer matrix. Section [3] translates the results of [TTU^ , ex- 
pressing spin form factors in terms of induced linear transformations, into the 
language of creation-annihilation operators. In Sectional we derive a number 
of identities relating trigonometric and elliptic parametrization of the Ising 
spectral curve. The main result is Lemma 14.21 which gives elliptic represen- 
tations for the matrix elements of A, B, C and D. Section [5] is devoted to the 
computation of two-particle form factors, i.e. the matrices D~^, BD~^ and 
D~^C. The principal results of this section are summarized in Theorem l5.6l 
Finally, the general factorized formula for the multiparticle form factors is 
established in Section |6l We conclude with a brief discussion of results and 
open questions. 



2 Transfer matrix diagonalization 

The two-dimensional Ising model on an Af x lattice is described by the 
following Hamiltonian: 

AI-l N-1 

-l3H[a] = ^ ^ (^ICxCrj,kO'j+i,k + l^v<^j,k(Jj,k+i^ ■ 

j=a k=o 

Spin variables {(Jj.k}j=o,....M-i take on the values ±1, and the boundary 

fc=o',...',Af-l 

conditions in each direction are either periodic or antiperiodic. This means 
that for j = 0, . . . , Af — 1 and fc = 0, . . . , — 1 we have cta/.S; = £xO'o,k and 
<^j,N = £yf j,o with Ex.y — ±1. To simplify some of the arguments below, it 
will be assumed that ICx.y G M>o. From physical point of view, this involves 
no loss of generality, as any real values of coupling constants can be made 
non- negative by a suitable change of spin variables and boundary conditions. 
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Typically one is interested in calculating n-point correlation functions 

{<^hM ■ ■■<^jn,k:} = ^^^^ ^-pH[a\ • (2.1) 

The 2^ X 2^ transfer matrix V can be thought of as an operator of discrete 
evolution in the x-direction. It acts in the space of states, which admits two 
equivalent descriptions: as a space of maps / : {1^2)^^ C or as an A''- fold 
tensor product (S) . . . (S) C^. To write V explicitly, it is common to define 
the operators 

Sj = 1 O (g) 1 (g) £7^ O 1 (g) . . . O 1, 

j times 

Cj = 1 (g) . . . (g) 1 (g) (Jj; (g 1 (g . . . (g 1, 

j times 

where j = 0, . . . , — 1 and Ox,v,z denote the Pauli matrices 
l\ /O -A fl 



The action of {sj}, {Cj} on the space of maps is given by 

i^jf) {^0^ ■ ■ • jCTiV-i) =crj/((To, . . . ,crjv-i), 
(Cj/) (o'O, • • • ,a-jv-i) =/(o-o, • • • ,crj_i, -(Tj,crj+i, . . .,crjv-i), 

where ctq, . . . , (Jn-\ € {1, —1}. The transfer matrix can now be written as 
y = (2 sinh 2/C^) ^ v} V^V} , (2.2) 

with 

. AT-l . . N-1 . 

V, = exp j /C* ^ Cj k Vy= exp j /Cj, ^ s^s^+i ^ . (2.3) 

Here Sjv = %so and /C* = arctanhe"^'^^ e M>o denotes the dual coupling. 

As already mentioned, V commutes with the global spin flip operator 
U = CqCi . . . Cjv-i- It also commutes with the operator T^^ of ^/-translations, 
defined by 

i'^Eyf) {(^0, ■ ■ • ,a-jv-i) = /(cTi, . . . ,aN-i,£yCro), 
as it may be easily checked that 

TeSjT-' = s,+u T.C.T-^ = Q+i, J = 0, . . . , iV - 1, 
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with Cn = Co- It is clear that \T^^, U] — 0, and therefore the operators T^^, 
U and V can be diagonahzed simuhaneously. Correlation function (|2.ip with 
Ji 5: j2 < ■ ■ • < Jri is equal to the trace ratio 



Tr 








Tr 







which can be straightforwardly computed in terms of matrix elements of sq 
in the corresponding basis of normalized common eigenstates. 

Remark 2.1 The transfer matrix considered in |llll26j is given, in our nota- 
tion, by 

N i i 

V = (2 sinh2/C:r)~ VyV^^ . 

Since V' is related to V in (|2.2I) by a similarity transformation, both defini- 
tions work equally well if one is interested in the partition function. However, 
Vx does not commute with the spin operators {sj}, and therefore V can not 
be used in the calculation of Ising correlation functions. This explains dis- 
crepancies reported in Section 9 of [TT] . 

On the other hand, the transfer matrix V' possesses the same set of eigen- 
vectors as the hamiltonian of quantum XY chain in a transverse field. Ising 
spin matrix elements between the eigenvectors of V should therefore coin- 
cide with the form factors of between the eigenstates of XY hamiltonian, 
which were found in [13] by the method of separation of variables. 



For j = 0,...,A^— 1, introduce the operators 



j times 



■ Cj" Sj , 



j times 



satisfying anticommutation relations for the generators of the Clifford alge- 
bra: 

{pj,Pk} = {q],qk} = 26jk, {pj,qk} = 0. 

It can be shown that 



V ^{2 sinh 2ICx 



1^ 



(2.4) 



where 



Af-l 



Vx = exp 



V = a,p, 
exp< —iJC 



and pn = —pa {pn = Pa) for — a (resp. v = p). Note that U anticommutes 
with all {pj}, {qj} and hence commutes with Va^p. 
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Consider discrete Fourier transforms 

N-l , N-1 

i 

Pe 



N-l N-1 

i=^e-^«,„ ,. = -1=5:^-^%. (2.5) 



J=0 ^ J=0 

Two sets of quasimomenta will be important for us: 0a = { ;^ 7 , ■ • ■ i 27r — } 
and 6p = {O, . . . ,2tt — In both cases one has inversion formulas 

P^ = ^T. * = ^ E (2-6) 

and anticommutation relations 

{Pe,Pe'} = {qe,qe'} = 259+6l\0 mod 27r, {^6,(79'} = 0. 

The conjugation by Va^p acts linearly on the generators {pj} and {qj}. In 
particular, Fourier components with 9 ^ 6^ {v ~ a,p) transform as 



VAZ\ 



cosh 79 ie "-^we \ ( Pe 
-ie^^w^e cosh 70 / \9e 

where the functions we and 70 are defined by 

we = cosh^ IC* sinh 2/Cy (l - ae''^) (l - /Se"^®) , (2.7) 
cosh 7e = cosh 2/C* cosh 2/Cy — sinh2/C* sinh2/Cy cos 0, 70 > 0, (2.8) 

with a — tanh IC* coth /Cj, and /3 = tanh /C* tanh ICy . It can be easily checked 
that 

sinh^ 70 = cosh^ /C^ sinh^ 2/Cy (l - ae'^) (l - ae''^) (l - ;3e*^) (l - /Je"*^) . 

From now on, unless otherwise is stated explicitly, we will work in the 
ferromagnetic region of parameters. The latter is defined by JC* < JCy or, 
equivalently, a < 1. In this case, it is convenient to write w±e in (12. 7p as 

w±e = 6±esinh7e, 

with 



- V (l-/3e»«)(l-«e-«)- ^^-^^ 

The branches of square roots in (|2.9p and below are defined so that they have 
a positive real part. 

Finally, introduce creation-annihilation operators of a- and p-fermions, 
which correspond to 9 Cz 9a and 6 Cz Op, respectively: 

hipl ^ e-'■'^^/b^P-e - i^/bleq-d, ^210) 
1 2^e = e'^^/bl^pe + i^/beqe- 
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Here —9 is identified with 2tt — 9 ioi 9 ^ 0. The operators in each of the two 
sets satisfy standard anticommutation relations 

Inversion formulas 

Pe = e-'^y6^(vle + V'e), qg ^ i^/b~g (yj^g - i^g^ (2.11) 

combined with Fourier transforms p.5p - (j2.6p . allow to express "antiperiodic" 
creation-annihilation operators as linear combinations of "periodic" ones and 
vice versa. 

These operators transform diagonally under conjugation by Va,p, 

Since {sj}, {Cj} (as well as {pj}, {qj}) generate full 2^ x 2^ matrix algebra, 
it follows from Schur's lemma that induced linear transformations fix Va^p 
up to a scalar multiple. Now taking into account that det Va^p — 1 and 
Tr Va^p > 0, we find that 

K = expi - ^ 7e (^eV'e " ^) f' = "'P- 

Fock states 

\9i,...,9k),^ijl...^ljvac),, 9i,...,9ke9,, (2.12) 
are therefore right eigenvectors of , with eigenvalues equal to 

I gee^ i=i J 

Here \vac)i, is a vector annihilated by all -06 with 9 G 0^. 

The operators {sj}, {Cj}, {pj}, {qj} are hermitian, hence p^ — p-g, 
q\ — q-g and the operator V'e is hermitian conjugate of ipg, just as the 
notation suggests. Introduce the dual vacua ^{vac\ = \vac)\, annihilated by 
all V'e with 9 ^ 6i,. Left eigenvectors of can then be written as 

i,{9k,. ■ ■ ,9i\ = \9i,...,9k)l = y{vac\^g^ ■■■■^g^, 9i,...,9k e 0^. 

We will impose the normalization condition i,{vac\vac),j = 1, which automat- 
ically implies orthonormality of the states (|2.12p . 

Recall that only a subset of eigenvectors of Va,p are actual eigenstates of 
V , cf. (j2.4p . For Ej, = 1 (i.e. periodic boundary conditions on spin variables) 
only the states with even number of particles are eigenvectors of the full 
transfer matrix. For Sy — —1, the number of particles in the surviving states 
of both types is odd. This follows e.g. from an explicit diagonalization of 
p in the limit /C* 0, and observation that the eigenvectors for /C* > 
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are related to those with /C* = by products of two-mode BogoUubov 
transformations, which commute with U. In fact since U also anticommutes 
with {il'g}, {"ipe}, it is sufficient to show that Z2-charges of the vacuum vectors 
\vac)a and \vac)p are equal to +1 and —1, respectively. 

To avoid confusion, in what follows we reserve the letters ipl, ipg for the 
creation-annihilation operators with 9 € 9a and denote the corresponding 
operators with 9 ^ 9p hy (fil, (pe- Our aim is to compute spin form factors 

J^m]n (^1 = a(^l, ■ • ■ , 0m\si\9[, . . . , 9'„) p = [p{9'„, . . . ,9[\si\9m, ■ • ■ , 0l)a] , 

(2.13) 

where 0i, . . . , 9m G Oa, 9'^, . . . , 9'^ e Op and * in the last expression denotes 
complex conjugation. As explained above, matrix elements (|2.13p are non- 
zero only if m and n are simultaneously even or odd. These are the only 
nontrivial form factors: since si anticommutes with U, its matrix elements 
between the eigenstates of V of the same type vanish. 

Remark 2.2 Let us show that the eigenstates of V constructed above diago- 
nalize the translation operator T^^ (since 79 = 7_e, the eigenvalues of V are 
degenerate and this is not automatic). Though translation invariance will be 
manifest in the final answer, which satisfies 

j-«„ (0, e') = e<^"- ^'-^"^1 ^0 j-w^ (0, e') , (2.14) 

it is instructive to explain it already at this stage. 

Indeed, the vacua \vac)a,p are non-degenerate eigenstates of the periodic 
Ising transfer matrix, and therefore they are also eigenvectors of T+ . It can 
be deduced from the limit /C* by continuity that the corresponding 
eigenvalues are both equal to 1. It is also straightforward to check that 



where e — ±1 and pn ^ epo, Qn = ^Qo- Performing discrete Fourier trans- 
forms, we obtain 

which in turn yields the desired formulas: 

T+\9u...,92k)a,p - e-^^Si"- |ei,.-.,e2fe)a,p, (2.17) 

T_ |0i , . . . , 92k+l)a,p = ^St' |0i , . . . , 92k+l)a,p- (2.18) 
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3 Spin form factors and induced rotations 

Just as in the above case of Va,p, conjugation by the spin operator s; acts 
Unearly on the Clifford algebra generators {pj}, {qj}- 

^^^^-f"-, J=0,...,N-1, (3.1) 

= sgn(/ -l~j)qj, 

where sgn(a;) = 1 if a; > and —1 if a; < 0. Recall that induced linear 
transformations fix si up to a scalar multiple. 

Let us combine {ipl}, {ipe} and {^Pg}, {'Pe} into iV-dimensional column 
vectors xp\ xp, tp'' , (p. Their entries are thus given by 2^ x 2^ matrices. It 
will be very convenient to write the result of conjugation of {y^J}, {'fig} by 
si as a linear combination of {ipg}, {V'e}' 

where A, B, C, D are N x N matrices. The entries of the 2iV-column in 
the l.h.s. satisfy canonical anticommutation relations, from which it can be 
deduced that 

AB^ + BA^ ^CD^ + DC^ = 0, AD^ + BC^ ^ 1. (3.3) 

It 

Since S; = Sj = , the operator s; is unitary. This yields further relations 

A = D., B = C. (3.4) 

Vacuum vectors \vac)i, {v = a,p) are fixed by the annihilation operators 
{ipe}, {^e} and the normalization ,y{vac\vac)i, = 1 up to inessential pure 
phase factors. This can be used to establish a relation between the two vacua. 



Lemma 3.1 We have 

\vac)p — A • sf^O\vac)a, (3-5) 

where 

|A| = |detZ?|^ = exp|-i ^ (^"'C), ^U. 

Proof. Lemma is formulated and will be proved under assumption that D 
is invertible. This point will be checked in Section [SJ where det D is shown 
to be non-zero by explicit computation. 

Acting on s'i^^O\vac)a with the A^-column y of annihilation operators, 
one finds 



<ps-^^0\vac)a = s^\Cip'' + Dip)0\vac)a = s^^OiCxp'' + DO-^ipO)\vac)a. 

(3.6) 
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For two matrices X, Y satisfying [X, [X,Y]] — one has e^Fe^^ = Y 
[X, Y]. From this and the fact that D^^C is skew-symmetric because of I 
we get 

O-^xpO ^ -D-^Cip' (3.7) 

Since iplvacja = 0, it follows from p.6p and p.7l) that (psY^O\vac)a — 0, 
and thus the vector s^^^Olvacja is proportional to \vac)p. 

Absolute value of the scalar multiple A in p.Sp can be determined from 
the normalization. Indeed, the unitarity of s; implies that 

1 = p{vac\vac)p ~ \X\'^ a{vac\0'' 0\vac) a- (3-8) 

To compute the vacuum expectation value on the right, one can e.g. write 
D~^C as G'^FG, where G is unitary and F has block-diagonal form, with 
2x2 skew-symmetric blocks containing nontrivial eigenvalues of D^^C and 
1x1 blocks containing zeros. Introducing new creation-annihilation operators 
= Gip\ ^ — Gxp which satisfy canonical anticommutation relations, and 
using that \vac)a is annihilated by ^, one obtains 

a{vac\0^0\vac)a = [det(l + FFt)]^ = 

= det (l + {D-^C) {D-^C^M ' = \deiD\-^. 

The last equality follows from the relation CC^ = 1 — DD\ which can 
be derived from Now ((XS)) implies that |A| = IdetD]^ and the 

lemma is proved. □ 

Since T+\vac)a,p = \vac)a,p, the factor A = a{vO'C\si\vac)p is independent 
of I and thus its phase can always be absorbed into the definition of vacua. 
We will therefore assume below that A € K>o. Under this convention, it 
follows from Lemma [3. II that 

a{vac\si\vac) p = IdetZ^I^. (3.9) 

Note that in the limit iV — oo the l.h.s. of this formula gives Ising sponta- 
neous magnetization. 

Let us now compute general form factors J-m^n{d, 6'). By (|2.13p and (|3.5p , 
one has 

= a{vac\i;g^ . . . Ve„ {Axp ^ + Bi(>)e'^ ...{Aip^ + Bi>)e'^0\vac)a, 

where we used (|3.2p to pass from the first to the second line. Now pull O 
in the last expression to the left vacuum, on which it acts as the identity 
operator. This can be done using p.7p and the fact that O commutes with 
all entries oi ip''. Moreover, (|3.3p implies that 
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where Z)"^ = {D^^)^, so that 

J-(i^„(0, e') = a{vac\{~D-^C'iP t + -0)^^ . . . {-D-^Ci, t + ^)^^ ^ 

Vacuum expectation value of any product of linear combinations of {ipl}, 
{ipg} can be easily computed using Wick theorem. The result is given by the 
pfaffian of the matrix of pairings between different combinations. In our case, 
there are three types of pairings: 

aivacli-D-'CxP^ +iP)e,{D-^'^^ + BxP)g,\vac)a ^ Dg\, , 
a{vac\iD'^',p^+B',p)g,{D-^tP^ + BiP)g, \vac)a^ {BD-') . 

Accordingly, we obtain 

Lemma 3.2 Form factors J^m.n{d,0') of Ising spin have the following rep- 
resentation in terms of induced rotations 



FliUe^e') = idetDi^ .pfi?, (^';' ^'j<'; ] , (3.10) 



x9 ^0' /.t 



where matrix elements of the blocks of skew- symmetric (m + n) x (m + n) 
matrix R are given by 



iRe>ce),k = {D-'C)^^g^, j,fc = l,...,m, (3.11) 
{Rexe')jk = -{Re'y.e)kj^ Dg^fi'^, j = l,...,m, fc = 1, . . . , n, (3.12) 

iRe'xe'),k = {BD-')g,g,^, j,fc = l,...,n. (3.13) 

Observe that matrix elements p.lll) - p.l3p coincide with normalized two- 
particle form factors: 

.(M-k.Kc), ^ (3.14) 
a{vac\si\vac)p ^ '^'^ 



a{vac\si\vac)p 
.{vac\s,\ee') (3.16) 
a{vac\si\vac)p ^ "^■'^ 

The reader might have noticed that the only important point for the 
above derivation is the formula p.2p . In particular, we have not as yet used 
the fact that s; belongs to the Clifford group, i.e. that ip\ ip and 
are related by a linear transformation. However, this does play a role in the 
proof of the following Lemma, which gives the explicit form of the induced 
rotation matrix. 
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Lemma 3.3 Matrix elements of A, B, C, D are given by 

„-t{i-k){s-e') 



Do 



2iNsin^ 



2iN sin ■ 



y%bg^ - 



(3.17) 
(3.18) 



where 9 G Op, 6' £ 6a. 



Proof. Express (fig in terms of pj, qj using (j2.10l) and the Fourier transform 
(j2.5p . compute the resuh of conjugation by si using p.ip . and rewrite it in 
terms of {ipl,}, {ipg'} using the inverse transform (j2.6|) and the formula (|2.1ip . 
Exphcitly, 



2si(pgS^ ^ = Si (^e'^^/b^pg + i\/bgqg^ S; ^ = 
= ^ H ^"'"^ i^^'^V^OSiPyS-^ + i^gSiqjS- 

= 7f E e"''' (e''y&I^sgn(? - j)p, + zv/6^sgn(Z - 1 - j)g,) 
]^ E E ^'''"''"^ (e^^ V^sgn(/ - j)pe' + *y^sgn(; - 1 - 

+ i^/b'gqg^ = 



2 ^ (.iiio'-e) 



2 ^ e 



Ibg, 



which gives C and I?. The formulas for A and B immediately follow from 
(lOll. □ 



Remark 3.4 An essentially equivalent result was obtained in fllj . Theorem 
4.1. However, the function bg appearing in the kernels of A, B, C, D differs 
from the one in [11]. This seems to be related to the inappropriate initial 
choice of the transfer matrix, see Remark 12.11 



4 Elliptic parametrization of the Ising spectral curve 

Let us now go back to the lattice dispersion relation (|2.8p . Introduce the 
variables z — e^^ , X = e^* and rewrite it as an algebraic curve 

A + A^""^ z + z^^ 
sinh2/Ca; [-sinh2/Cj, = cosh 2/Ca; cosh 2/Cy. (4-1) 
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Topologically this is a torus realized as a two-fold covering of P^, branched 
at 4 points 

zi,2 — (tanh/C* coth/Cj,)^"'^ ~ a^^, 
Z3,4 = (tanh /C* tanh/Cj^)*"'^ = (3^^. 

The points /3 < a < oT^ < /3^^ can be mapped to < fc < k^^ < cx) by a 
Mobius transformation. Such transformations preserve the anharmonic ratio 

(zi, Z2; Z3, Z4) — -, and therefore 

[Z2 - Z3){Zl - Z4) 

13 -a sinh 2/C! 

fc = 



al3-l sinh 2/C 



y 



The parameter fc plays the role of the modulus of the Jacobi elliptic 
functions uniformizing (|4.1I) . One also needs to introduce elliptic integrals 



^il-t^)il-kH^)' Jo ^il-t^){l~k'H^) 

where the complementary modulus is defined by fc'^ = 1 — fc^. It is known 
(see Theorem 2.2.1 in [3S]) that the functions 

ziu) = ^^^j^, (4.2) 
sn [u — IT]) 

X{u) — [k sn {u + irj) sn {u — irj)] ^, (4-3) 

with u G [—K, K) X i[—K\ K'), provide a uniformization of the Ising spectral 
curve (|4.1I) . The real parameter rj G ' determined by sinh 2ICx = 

i sn 2ir] (it is related to a of j25j by a = 77 + ■^). 

Elliptic parametrization (j4.2p - (|4.3p establishes a one-to-one correspon- 
dence between the "physical" cycle Ce = {{z, A) = (e*^ e'^«) | 6* e [0, 2n)} and 
the set C„ = {m I Rew e [—K, K), Imu ~ 0}. It is convenient to denote by ug 
the image of the point {e^^,e''^) G Cg in C„. In particular, one has uo = —K, 
Wtt = and, more generally, ug = —U2TT-g for G (0,7r]. 

The rest of this section is devoted to the derivation of a suitable elliptic 
representation for the elements of the induced rotation matrix. We proceed 
by establishing a number of identities between the two parametrizations. 
First note that for 9 G [0, 2n) holds 

e 2 = -Vk sn{ug ^ irj). (4.4) 



Indeed, the squares of both sides of this relation coincide because of 
(|4.3p . Since they never vanish, it is sufficient to fix the sign for one value of 
9. Set 9 — n, then the l.h.s. becomes ^pze"^"/^ and sn{ug =F irj) in the r.h.s. 
reduces to ^psniry. On the other hand isnii] > 0, and the result follows. 
Setting u = Ug, u' = ugi , v ^ irj in the easily verified formula 

sn(u -I- v) sn(u' — v) — sn(u — v) sn(u' + v) , /n „ 

— -sn(u - u ) sn2u. 



1 — k'^sniu -\- v) sn(ii — v) sn(u' + i')sn(u' — v) 
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and using (|4.4I) . one finds that 



sn{u0 -ue')^ sinli 2/C,y 7,+7,, • (4-5) 

This identity has several important consequences. Notice especially that for 
9' — 0,n,2Tr — 9 it reduces to 



sn(ue + K) sin ^ 
sinh2/Cj, ^sinh^:^' 

sn ug cos I 



sinh2/Cy sinh2ii±2i' 

sn 2ue sin 9 

sinh 2/Cj, sinh 7e 



(4.6) 
(4.7) 
(4.8) 



Since cosh 79 — cosh 79' = sinh2/C* sinh 2/Cj, (cos 6*' — cos 6*) by (|2.8p . it can be 
inferred from (|4.7p that 

sinh ^ ^ 

sinli T sinli '^^ '^'^o' 

1 - fc2 sn^ug sn^ug, = ^- — . (4.10) 

sinh sinh ^ ' 

Next we want to compute cnug and dnug. For that, set 9' — 9 in the 
last identity and combine the result with the doubling formula sn2u = 
2 sn u cn u dn u 

-T. — -. and (14.7P to fix the product cnuednMg. Together with (14. 6p 

1 — fc^sn^u 

cn u 

(since sn(M + K) = and cnug, Anug > 0), this gives 

dnu 



. 9 / sinh 2/C„ sinh 7^ / sinh 7^ sinh '^0+^'" 

cn Ug — sm — 4 / 1 — ] , dn Ufl = 4 / . 

2 Y sinh 2^ sinh 2^ V sinh 2/Cy sinh 

(4.11) 

To find an elliptic representation for bg in (j2.9p . recall that 

be = cosh^ JC* sinh 2/C„ ^. 

" smh7e 

As a function of 9, the numerator of this formula is a linear combination of 
1, COS0 and sin 9. On the other hand, (j4.9p and (|4.10p imply that 

2fc sn^ue _ cosh 7^ - cosh 79 (4 12) 



1 — fc-^sn'*ue sinh 7.n- sinh 7e ' 
1 + fc^sn'*ue cosh 7^ cosh 79 — 1 
1 — fc^sn^ue sinh 7^ sinh jg 



(4.13) 
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Hence bg is a linear combination of the l.h.s.'s of (j4.12L (|4.13p and (j4.8p . A 
little calculation then gives 



be = ik sn 2ug 



2fc sn ue I + k sn ug (dnug + iksnug cnug) 



1 — k'^sn^ug 1 — k-^sn^ug 1 — k^sn^ug 

Taking the square root and fixing the signs via ^/b^ — 1, one obtains 
/— 1 dnug ±iksiiugciiug 



(4.14) 



From (j4.14p we can now deduce what we really want — namely, a nice 
elliptic representation of the matrix elements of A, B, C , D in Lemma 13.31 



Lemma 4.1 Let 9, 9' e [0, 27r) with 9±d' ^0, 2n. Then 



sinh 



dn(w0 - ug,) 



2sin^ 
1 



2sin«±^ 



T Vbgbg' - 



1 



^/suiK^yg'suiil^ Sn{ug — Ugi) 

—i sinh 2/C 



(4.15) 



Vbgbg^ J y/ siuh Jg siuh 79' 

Proof. From (j4.14p it follows that e.g. 



cn{ug - ug,). (4.16) 



1 



2 sin 



Standard addition formula dn(M — u') 



dn ug dn ug> + fc^sn cn ug sn ug' cn ugi 
sin »^ ^(1 - fc^sn^ue) (1 - k^sn^ug,) ' 

dn u dn u' + /c^sn u cn m sn u' cn u' 



1 — fc^sn^u sn^? 



transforms the right side of the last relation into 

1 — k'^sn'^ugsn^ugi dii{ug — ugi) 

V'(l - Psn^ue) (1 - k^sn-^ug,) sin ^ 

sinh 2i±Iei 

By (|4.10p . the first factor may be reduced to 



-y/ sinh 7e sinh jgi 

sired identity (|4.15p becomes an immediate consequence of (|4.5 
Similarly, applying (j4.14l) to the l.h.s. of (I4.16P we get 



and the de- 



2sin^ 



7 \/bgbg, 



ik 



sn Ug cn Ug dn ugi + dn ug sn wg' cn wg' 



^/hbi' J v'(l - fc^sn^Me) (1 - A;2sn4„e') ' 



Now (j4.16p can be obtained from the addition formula 

, , snucnudnu' + dnusnu'cnu' 

cn(M — u ) sn(u + u) = — 

in a way analogous to the above. 



1 — k'^sa^u SB?i 



□ 
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Finally, let us introduce two diagonal matrices (y = a, p) with elements 

(^,),,, - ^==5e,e', 0,6' ee,, (4.17) 
Vsmn 76/ 

and two non-diagonal matrices <P, ^ defined by 

^ dnK^^^ tf.,,, ^cnK-ueO, eeep,e'eea. (4.i8) 

Combining the results of Lemmas 13.31 and 14.11 we arrive at 

Lemma 4.2 Matrices B, C and D defined by \3.11^ - f3.1^) satisfy the fol- 
lowing identities: 

—iN 

Idet D\ = , (sinh2C.)"|d.t^| 

ylleeep sinh7e Heee^ sinh7e 

Proof. Check that B^C^ J.inh21Cl ^ ^ _sinh2^ A^^AaH 

N ^ iN ^ 

Thus in order to obtain all finite-lattice form factors of Ising spin, it is 
sufficient to compute the inverse and two matrix products 'F<P^^, 'P~^'P. 
This will be done in the next section using that <P, 'I' defined by (|4.18p are 
special cases of elliptic Cauchy matrices. 



5 Elliptic Cauchy matrices and two-particle form factors 

5.1 Elliptic determinants 

Consider 2N complex variables xi, . . . , z^v, ■ • ■ , Vn- The determinant of 
the Cauchy matrix Vij — {i,j — 1, . . . ,N) is a rational function of 

{xi} and {yi}, which vanishes whenever Xi = Xj or yi — yj for some pair of 
distinct indices i,j. It is well-known that there are no other zeros and the 
determinant is given by 

There is an elliptic version of this relation due to Frobenius [8] . Namely, if 
we denote by 'di,,,i{z) Jacobi theta functions of nome q = e'^*'^, then for the 
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-diixj - Vj + a) 
di{xi - yj)-di{a) 
parameter a one has the identity 



eUiptic Cauchy matrix Vij = — tTw — n" depending on an arbitrary 



^1 (Ef - Ef y^ + «) nf<, ^1 - ^3 ¥i ivj-y^) ,,,, 

detV ^ ^ ir-r-^ iv • (5-1) 



Corollary 5.1 We have 

Vmn = X (5.2) 

'^l (Sf Xi - Vi + "j '^l {^n - Vm) 

^ llf ^1 {Xn - Vi) nf ^1 - Xj) 

Vu^n ^1 i^n - X{) n,^„ (y™ - Vi) 

Proof. Use (|5.ip to compute the determinant and cofactors and (|5.2p fol- 
lows. □ 

Theta functions are related to the Jacobi elliptic functions by 
snu = , „ , , cn M = , „ , , dn M — 



(5.3) 

where -di = i?i(0) for i = 2,3,4. Elliptic modulus and half-periods are given 

by 



q 

n 

The functions t?2,3,4(2) can be obtained from ^i{z) by shifting its argument, 



/fc = -I, 2K = TT-dl, 2iK' = 7rr?9?. 



^3(^)=e-(^-^)^i(z + |-^), (5.5) 
i94(^) = ie~<"-T)^^ " t) ■ ^^'^^ 



Therefore one can deduce from (|5.ip . (|5.2p the determinants and inverses 
for a number of matrices with entries written in terms of the Jacobi elliptic 
functions. 

Lemma 5.2 For 1 < i, j < N with even N we have 

/ N \ 2 

det (^Vk STi{ui — Uj)^ — I JJ^ Vk sn(wi — Uj) 1 . (5-7) 

\i<j ) 
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Proof. Set Xi — -d^ "^Ui, iji — -d^ '^Ui + ^ and use (|5.3p - (l5.6l) and the formula 

(^) = ie-T^i?4 to write 

The determinant on the left side of (|5.7p can therefore be deduced from (|5.ip 
with a = Since in our case yi = Xi + ^ and in particular J^f -^^ ~ 
Vi = it reduces to 

iNnT 2^'"'^ ) TT- di{xi - Xj)'di{Xj - Xi) 



i<j 

To transform the thcta functions back to sn's, use again (j5.3l) - (j5.6l) and the 
quasiperiodicity relation ■di{z + I-kt) = (— l)^e^*'^^ "^"^^^^i^i (z). which holds 
for ^ e Z. ' □ 

Remark 5.3 Taking the square root at both sides of (|5.7p and fixing the signs 
via the rational limit, one obtains a pfaffian version of this formula: 

N 

Pf ^VA: sn(Mi — Uj)^ = JJ^ Vk sii{ui — uj). (5.8) 

We note that this identity has already appeared in the proof of Theorem 5.0 
in [271, where it was used to obtain multiparticle Ising form factors on the 
infinite lattice from the two-particle ones. 

Lemma 5.4 Let denote a matrix with elements <Pij = — - — ^ with 

sn[Ui — Vj) 

i, j ~ 1, . . . , N . Then 

J ■&2'^4 _iCy-«^._y-« xr\ /sr^N „W tt ttt\ 



Uf<j "^li^i - Xj)^iiyj - Vi) ,^ 

X ^ , (5.9) 



mn 



where Xi = 'D^'^Ui, yi = 'D^'^Vi. 
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Proof. From (j5.3l) and the relation t^i (f — ^) = e "•^ "dz it follows that 
^i," can be written as 



Now to derive ([521), (|5?T0| . it suffices to set a = f - ^ in ([5?I]), dO]). □ 



5.2 Theta functional interpolation 

Let zi, . . . , za/, /i, . . . , /if be 2M complex variables. The Lagrange polyno- 
mial P{z) is the unique polynomial of degree < Af — 1, interpolating the 
set of points {zi, fi) (in other words, fi = P{zi) for i — l,...,Af). As is 
well-known, 

M M 

... -^2 ^1 

Introduce M additional parameters and the (Af — l)th degree 

polynomial 

M M 

f{z) = l[iz-z',)-l[iz~z,). 

i i 

If we now set fj = /(z^) for j = 1, . . . , M, then the evaluation of the coeffi- 
cient of z^~^ at both sides of (|5.1ip leads to the identity 



_ M M 



The last relation has an elliptic analog, which turns out to be very helpful 
in the calculation of the products BD"^ and D'^C . Namely, if the parame- 

M M 

ters {zi\, {Zj'} satisfy the balancing condition X] -^i ~ S -^i — 0: then 



Y ^ —=Q. 5.12 

This identity appears in a somewhat disguised form in [35], Example 3 on 
p. 451, where it is formulated in terms of the Weierstrass cr- function. It can 
be proved by induction on M , using the Riemann's addition formula in the 
base case, see e.g. Theorem 7 in [3T]. Another proof, based on the Frobe- 
nius determinant formula, is given in the Appendix of [12]. The identities of 
type (j5.12l) arise in [211122] in the analysis of a discrete time version of the 
Ruijsenaars-Schneider model. 
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Lemma 5.5 Let ^ denote a matrix with elements — cn{ui — Vj) (i,j 
1, . . . ,N) and {xi}, {yi} be as in Lemma\5.4\ Then lue have 



11 11 ■ ^'-''^ 

Proof. We will prove only the first identity, since for the second the argu- 
ment is completely analogous. By (j5.3p . one can write 



^2 ^1 {Vm -Xi-^] 



Lemma [53] then implies that the l.h.s. of (j5.13D is given by 



^1 (Ef - Ef + f - n.",. - -0 



N 



(j/m " + f ) n»^n ^l(?/m " Xj) ' 

^1 n;^™^i(y™-2/.). 



If we denote 



2 ' 



Zi = yiioi i = l,...,N, z^+i ^ xi + _ 
z[ = Xi for i = 1, . . . , n — 1 and for i = ri + 1, . . . , A^, 

•^n = -^n ~ Ei •^i ~^ Ei J/i + § + ^ ^-nd ^jv+l ~ ~ i 

then the sum in (15.151) can be written as > — ^rrr^ . This almost 

coincides with the sum in the l.h.s. of (j5.12l) with M = + 1. The missing 
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{N + l)th term is given by 



2 . 



Since the balancing condition is satisfied, our sum is equal to minus this 
missing term. Together with the first line in (jS.lSp . this leads to (|5.13p . □ 



5.3 Trigonometric reduction 

The matrices ^ and corresponding to the Ising case, are indexed by two 
sets of quasimomenta. Explicitly, Uj = ug with 9 = ^^''^ S 9p, and Vj = ugi 
with = {i ~ \) G 0a- Since for G (0, 27r) each set contains both 9 and 
27r — 0, there are additional constraints on {xj}, {yj}. Their form slightly 
differs for odd and even values of N: 

— For odd iV, the value 9 = belongs to the periodic set of momenta, 
but = TT is in the antiperiodic one. Then for all j = 2, . . . , and 
k = 1 , . . . , -^^^ we have 

TT 

Xj + XN+2-J ^yk+ VN+l-k = Xi + - = y{N+l)/2 = 0. (5.16) 

— For even N, both exceptional values 9 = 0,tt belong to the periodic set. 
Thus for all j = 2, . . . , Y and fc = 1, . . . , one has 

TT 

Xj + XN+2-j =yk + yN+l-k = Xi + - = Xn/2+1 = 0. (5.17) 

In particular, in both cases we have 

N N 

Y^x,-J2y, = -^. (5.18) 

i i 

This allows to rewrite (|CT)) - (|ETU)) and ^^J^-^J^ as 

fdet^>)^ = -^^^-^ <p-^ = (5 19) 

^ ^ nf/.5.' ™ sn(.„-.„)' ^^-^^^ 



{'F'P ^)^^=fnh{xi + —)an{ui-un), 



h iyi 



—-sn{v,n ~Vl), 



(5.20) 
(5.21) 
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where h{z) = TT — ^ and {/„}, {gm} are defined by 



Notice tliat (jS.lSp . combined with the quasiperiodicity of "ffiiz), also im- 
plies that h{z + ttt) = —h{z). Further, it is easy to check using (j5.3p . (I5.6p 
that the quantities Xn = —fnh{xn — 3-) and Km = gm/h (^i/m — ^) are 
given by 

Thus, to write (|5.19p ~ (|5.2ip in terms of sn's, it sufhces to find a suitable 
representation for the function ^{z,z') = h{z)/h(z'). For that, one can use 
the constraints ([5l^ - (ICT7)) . 

Suppose e.g. that N is even, then 

(5.23) 

where we have put together the theta functions with arguments coming in 
pairs according to (j5.17l) . In each product of two theta functions, use the 
well-known addition formula, written in the form 

After substitution into (|5.23l) , almost all of the theta functions and denomi- 
nators appearing in (j5.24p cancel each other, so that fi{z, z') reduces to 

cniJgZ (1 -I- sniJgz') JX- (sn??!^ — snui) (snjJgz' — snwi) 
cndlz'{l + sndlz)\}^ {sn^lz - syiv,) {siidlz' - snu,) ' ^^''^^^ 

Similar manipulations for odd N lead to the same formula. 

If we now denote A(m, w) = /i {-d^^u - ^, i^g - then by 

, , dnu (1 -I- fc snw) -1^ (1 — fc snuj sum) (1 — fc snwi snw) /r 

Xiu.v) = ; r ; -— ; r . (5.26) 

dnu (1 -I- fc sum) -"-^ (1 — K snvi snu) (1 — k snu^ sn-y) 

Relations (|5.22p and (|5.26l) then allow to rewrite (|5.19p - (|5.2ip in terms of 
the Jacobi elliptic functions: 



i^-^^f- y,N^' \; ^ y ^- = '3:t;7 - (5-27) 



(-1)^\/1 - _ KmXnK'"m,Un) 
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(tf'^> =x„A(u;,w„)sn(u; -u„), (5.28) 

{^-^^)^^^ = KmX{Vm, Vl) Sn{vi - Vm)- (5.29) 

So far we have used only a reflection symmetry of the periodic and an- 
tiperiodic spectrum of momenta. That these momenta take rather special 
values (recall that e'^ is an iVth root of ±1) leads to further simplifications. 
It will be convenient to switch the notation and write all matrix indices in 
terms of momenta, as suggested in the beginning of this subsection. Thus, 
e.g. 



lle'eep,e'5^esn(ue - US') 1 le'ee„,e'#e " 



where in the first formula 6 ^ Op and in the second one 9 ^ 0a- Likewise, for 
e, 9' e [0, 27r) denote 

, . dnug T-r 1 — ksaugii snug -p-i- 1 — ksnug// snugi 
\e,e' = \[Ue,U8') = II -; ; II ; . 

dnufl/ 1 — fc sn Mo" sn Mfl/ 1 — Ksnufl"SnMfl 

0"eep o"eea 

(5.31) 

Our aim is to express xe, and Xg,g' in the initial trigonometric para- 
metrization (i.e. in terms of 9 and 7e). Introduce the function 



Ylo'aa sinh 24221 

i.g=\nt±l^ ;r4^, 6l€[0,27r). (5.32) 

TT sinh ^0+^0' l : ; \ j 



Use (|4.1ip to rewrite the first factor in (I5.3ip and transform the remaining 
two products with the help of (|4.10p . The calculation is slightly different for 
odd and even N, but the final result is given by the same simple formula: 

Ae,e' = e(''«'~''«)/2^ (5.33) 
Similarly, the identity (j4.5p allows to write (j5.30p as 



_,^ sinh2/C, rie.ge^sin^ f n^n o,^ 
Xe=e "^-r . g-O' for 6* e 0p, 5.34 

sinh2Ar TTfl'co sin 4^ 

^,=e-^ smli2/e for^?e0,. (5.35) 

smh79 Yl^,^^^^^,^^sin!^ 



25 



The products in the r.h.s. can be easily computed exphcitly. One has 



2^-^ n sin ^= (-1)^-1 sin ^, (5.36) 
e'eOp 

2^-1 H sin^= (-1)^ cos^, (5.37) 



2^-1 H sin^=(-l)^-iiVcos— ioveeOp, (5.38) 



2 ' ' 2 



2^-1 n sin^=(-l)^-iiVsin^ ior 6 e 6^. (5.39) 



The first two relations are valid for any 9 £ C and follow from the obvious 
product formulas Yie'eOp^"^ ~ e"^^ ) — — 1 and Yie'eOa^'^ ~ e^^ ) — z'^ + 1. 
The third and fourth formula may be deduced from (j5.36l) - (j5.37p by taking 
appropriate limits. 

Substituting (|5.36p - (l5.39p into (I5.34p - (|5.35p . and combining the result 
with (|i?5|) and we can rewrite ((07)) - (IOg)) as follows: 



Af2Af /-i _ 1.2 

(detg>)^= ^ n e''^/'sinh7e H ^-^^^^ sinh7e, (5.40) 

(smh2/Cj,) g^g^ 



sinh2A:„e(''''-'^<'')/2 sinh 2:^ 



-7s' 



iV2sinh7esinh7e' sin ' ^e0a,f?'e0p, 

(5.41) 

sinh^ 2Ar sin 
V hfi' iVsinh7e sinh^^' ' - ^ > 



Now from (j5.40p ~ (l5.43p and Lemma W?^ one easily derives the main result 
of this paper: 
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Theorem 5.6 Let B, C and D be the matrices defined by {3.11\ l- ^.18\) . 

Then 



BD- 



|detD| 



sinh^^ 



N^J sinh 7e sinh 79/ sin 



(5.44) 



sm ■ 



sinh 2JCx N^sinh-ie sinh 79- sinh 2^±2ai 



(5.45) 



-Hi-h)(e+e') sinh2/C,; 



,(i/e+^e')/2 



sinh 2/C2; 7VVsinh7e sinh 79/ sinh 2»±2£i 



(1 - fc^) n n 



(5.46) 
(5.47) 



w/iere 6* e Oa, 9' G 6»p m ^5^W^, 9' e 6/p m [5^I\ ), 9, 9' £ Oa in and 
the function vg is given bv S5.32\) . 

By p.9p and p.l4p - p.l6l) . Theorem 15.61 gives the vacuum expectation 
value a{vac\si\vac)p and two-particle form factors of Ising spin. Multiparticle 
form factors are determined by Lemma 13.21 



6 Multiparticle form factors 



We now compute general matrix elements Tm.,n{d,0') using a trick learnt 
from 11 . There it was shown that, roughly speaking, if the conjectures of [SI 
[7j are true for two-particle form factors, then they also hold for multiparticle 
ones. Theorem 15.61 implies that the last statement is no longer conditional. 

The argument of [11) may be sketched as follows. Introduce an (to -t- n) x 
(to + n) diagonal matrix Q with non-zero elements defined by 



' exp{-z(;-l)g,+^.e72} 
sinh 7ej 

cxp{z(/ - |)g;._„, - ;.e.._^y2} 

y/iV sinh 7e^_^^ 



for 7 = 1, . . . , TO, 

for i = TO + 1, . . . , TO -|- n. 



Taking into account Theorem 15. 6[ the skew-symmetric matrix R in Lem- 
ma 13.21 can be written as 



R = -ipQRfl, 



R 



Rexe Rexe' 
Re'xe Re'xB' 
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where p = \ —. — ^ and matrix elements of the blocks of R are given by 
smh2/C:r 

(^«xe)jfc = — J,k = l,...,m, (6.1) 



sinh 
sinh 



2 

78 ,-+79 



p sin ^ ^ '° 



■2) 



(7 ■ — [/ ,, 

/? sm ' ^ 



(^^'x.'),. - . > j-/c = l,...,n. (6.3) 



sinh 



12 Ifc. 



Let us recall the identity (|4.5p from Section 21 which we rewrite in the 
form 

Vksn{ug-Ug,) = [Vksn{ug-Ug,±lK')y^ = ^' ^' ^ t^' 

(6.4) 

Introducing (m + n) variables 

forj = l,...,m, 
I wg' tor j ~ TO + 1, . . . , TO + n, 

matrix elements of R can be written as Rij — \/ksn{ui — uj) with = 
1, . . . ,TO + n. The pfaffian (I3.10p can then be easily computed using the 
elliptic pfaffian identity (jS.Sp : 

m+n m+n 

ViR = i-ip)"^ -detn •Pfii= i-ip)"^ Yl ^" 11 ^sn{ut~uj). 

i—l i<j 

Transforming sn's back to trigonometric functions by (j6.4p - (j6.5p and taking 
into account p.lOp and (|5.47p . we finally obtain the desired general formula 
for finite-lattice form factors of Ising spin: 

Theorem 6.1 Spin matrix elements J^m,n {0,0'), defined by \2.13\) . for even 
m + n can he explicitly written as 

2mn— "^i'" 



/sinh2X;„\^ g-.(;-i)9,+.«,/2 JL ^(i-m-ed2 



7ej +7e'. 



\ sinh 2K,x 

sin T-r sin -n- sinh 



ii . -,78, +78, ii 7e'+7e' ii 



l<i<j<m 2 l<i<j<n Smh — ^-^5 — - l<i<m Sm ■ 

i<i<" 
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where ^ = |l — (sinh sinh 2/Cj,) 
and 



the function vg is defined by (5^ 



n 



'Ug/4 ^ 



sinh 



2 78+78' 



'ee, 



sinh 



79+7a' 
2 



'ee, 



sinh 



78+70' 
2 



(6.7) 



This resuh coincides with the conjectures of [Bl[7]. It clearly satisfies trans- 
lation invariance constraint (I2.14p . In the thermodynamic limit — > 1, 

^ Q, the spectra of quasimomenta become continuous and (j6.6p repro- 
duces infinite-lattice form factors found in [271129) . In particular, considering 
the N 00 limit of (|5.47p . it is straightforward to recover Yang's formula 
("') = for the spontaneous magnetization |33) . 



Remark 6.2 We finally comment on the paramagnetic region of parameters 
(/C* > ICy). The above results remain valid if one performs analytic con- 
tinuation in e.g. /C*. Nontrivial dependence of (|6.6p on /C* is hidden in the 
functions 79 with 9 G 9a,p- The result of analytic continuation of almost all 
such 7e from the region /C* < ICy to /C* > ICy coincides with the definition 

-70. Then one finds 



\. The only exception is 70, which is continued to 
that under continuation 



sin ■ 



sinh 2:4^ 



sinh2/C:, sinh^i^ 
sinh 2/C„ sin % 



p!^e/2 



/sinh2/C:, sinh^i^ 
sinh 2/C„ sin ^ 



for e e (0, 27r) and 





/ sinh 2JCx 


( e-''"/^ ^^ 


sinh 7o 


Y sinh2A:j, 


^ ^JN sinh 70 / 



smh2K.y 



sinh 2JCx N sinh 70 



(6.9) 
(6.10) 



Quasiparticle interpretation of eigenvectors also changes. If the particle with 
9 = was initially present or absent in a p-eigenstate, then it is annihilated 
(resp. created) after analytic continuation. This means e.g. that the number 
of particles in p-eigenstates of the periodic Ising transfer matrix becomes odd 
instead of even. Relations (|6.8P " (I6.10I) then show that spin form factors are 
given by the same formula (|6.6p . with odd m + n and appropriately adjusted 
first numerical factor. 

It is also possible to derive this result along the above lines. However, for 
/C* > ICy the definition of the creation-annihilation operators (|2.10p should 
be modified and the corresponding matrix D is degenerate, which leads to 
additional subtleties. 
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7 Discussion 

In this paper, exact expressions for finite-lattice form factors of the spin 
operator in the two-dimensional Ising model were obtained. Starting point of 
our derivation was the idea of expressing induced linear transformations of 
fermions in a particular basis, which was put forward by Hystad and Palmer 
in [Tn[26] . New crucial ingredient is the use of elliptic determinants and theta 
functional interpolation. The presented approach seems to be quite natural 
and straightforward. It is likely to be applicable to other free-fermion lattice 
models, such as those considered in P^ITSII^ . 

A more ambitious challenge is to complete the program of form factor 
derivation for Zjv-symmetric superintegrable chiral Potts quantum chain, 
which presumably possesses a hidden fermion structure. Two parts of the 
hamiltonian of this model generate Onsager algebra and the space of states 
decomposes into a set of invariant irreducible subspaces (Onsager sectors) 
with Ising-like hamiltonian spectrum in each of them. 

Using Baxter's idea of extending Onsager algebra by the spin operator 
[2], its form factors between the hamiltonian eigenstates were found in [T3] 
up to unknown scalar factors for each pair of Onsager sectors. In the case 
N — 2 (quantum Ising chain in a transverse field), Onsager algebra can be 
embedded into a fermion algebra, so that all irreducible representations of 
the former are combined into one irreducible representation of the latter. 
This allows to fix all unknown scalar factors and obtain form factors of the 
quantum Ising chain |12j . It is therefore natural to look for a fermion-like 
algebra extending Onsager algebra for general N. 
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